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Abstract. We consider a fermionic dark matter candidate interacting via a scalar mediator 
coupled with the Standard Model through a Higgs portal. We consider general setting includ¬ 
ing both scalar and pseudoscalar interactions between the scalar and fermion, and illustrate 
the relevant features for dark matter abundance, direct search limits and collider constraints. 
The case where dark matter has a self-interaction strength (cry) /m^ ~ 0.1 — Icm^/g is 
strongly constrained, in particular by the Big Bang Nucleosynthesis. We show that these 
constraints can be alleviated by introducing a new light sterile neutrino N. The allowed 
region for the extended model consists of a triangle at lO(sin0)“^/® MeV < rriN < 1 GeV. 
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1 Introduction 

Recently there has been a lot of interest in models with extended matter sectors very weakly 
coupled with the Standard Model (SM) fields, see e.g. [1-11]. These are mainly motivated 
by the need to explain the dominating dark matter (DM) component of the energy density 
of the Universe and simultaneously hiding the beyond SM fields from the direct searches at 
collider experiments. Alternatively, the extra matter fields coupled with the SM scalar sector 
also affect the properties of the electroweak phase transition and enhance the possibility of 
successful electroweak baryogenesis. 

The simplest possibility of a hidden sector consists of a single singlet scalar coupled 
with the SM via the Higgs portal, see e.g. [1, 12-14] for DM studies and [15, 16] for aspects 
of electroweak phase transition. While the singlet scalar can in principle act as DM and lead 
to strong first order phase transition, in light of present constraints for the DM abundance 
and from the direct search experiments, these features cannot be realized simultaneously 
[14]. To alleviate this issue one is led to consider more complex hidden sector structures. 
The possibilities of singlet fermions and vectors have been considered in the literature [5-11]. 
For example in the case of singlet fermion DM interacting with the SM via Higgs portal one 
can easily provide sufficient DM abundance and simultaneously use the scalar mediator to 
provide for strongly first order electroweak phase transition [8]. 

Further motivations for extended hidden sectors arise from the possibility of DM self¬ 
interactions as has been studied in various contexts [17-37]. Such interactions would allow 
to reconcile the tensions between the observations and simulations of small-scale structure of 
collisionless cold DM (CCDM): First, as the observations of dwarf galaxies confirm, the DM 
density profile is not cusped towards the center of the galaxy as seen in simulations of CCDM 
[38] but exhibits a flat core [39, 40]. Second, the so-called missing satellite problem [41-43] 
and the too big to fail problem [44] can also be explained by sizable interactions between the 
DM particles, although alternative astrophysical explanations also exist [45-51]. 
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In this paper we consider a simple paradigm for self-interacting DM: a singlet fermion 
ijj interacting via scalar mediator s. We consider in particular the parameter region where 
the scalar mediator s is light. In [22] it has been pointed out that self-interactions of a 
DM particle mediated with a light force carrier s and satisfying (m^/10 GeV)(ms/100 
MeV)^ ~ 1, can provide the flat density profile around the core of dwarf galaxies for a wide 
range of interaction strengths between ijj and s while evading constraints on self-interactions 
from galactic and cluster scales. However, there exist additional stringent constraints. Most 
notably, if the scalar s couples to the SM via the Higgs portal [2, 52-58], it must decay before 
the onset of the big bang nucleosynthesis (BBN) around t ~ 1 sec. in order that the success 
of the BBN explaining the abundances of light elements is not endangered. This sets a lower 
bound on the strength of coupling between the hidden and visible sectors which is too strong 
[31] in light of present constraints for the invisible Higgs decay at LHC and direct search 
constraints from the LUX experiment. Generic ways to evade this problem were proposed in 
[59]. 

Our paper can be viewed as a significant complement to earlier literature on singlet 
fermion DM [5-8, 59]: We investigate the effects of both scalar and pseudoscalar couplings 
between the DM fermion and the scalar mediator. Via Monte Carlo scan, we probe wide mass 
scales for the singlet sector especially focusing on the range where the DM self-interaction 
cross section is large enough to reconcile the observations and simulations of small scale 
structure. In addition, we take systematically into account the phenomenological constraints 
from DM abundance and direct detection as well as theoretical constraints of stability and 
perturbativity. Furthermore, we confirm that the BBN constraint can be resolved in this 
model via the mechanism proposed in [59]. 

The paper is organized as follows: In section 2 we define the model and outline the basic 
computations constraints on DM abundance and direct detection. In section 3 we study the 
strength of DM self-interactions and its relation to the mass of the scalar mediator. The 
consequences of the existence of a light mediator on the physics of the early Universe is 
discussed in section 4. In section 5 we present our conclusions and outlook. 


2 Model 

2.1 Scalar sector and the dark matter Lagrangian 


We consider an extended scalar sector described by the potential 

2^ ' 3 


^(</’) s) — 




( 2 . 1 ) 


which provides typical scalar portals between the hidden singlet sector and the SM Higgs. 
The field (j) is the usual SM Higgs doublet and s is a real singlet scalar. The Higgs doublet 
is written in terms of its components as 


^ (73(^ + ^ + w))’ 

where v is the vacuum expectation value which at T = 0 has the value v = 246 GeV. The 
stability of the tree level potential requires that 


^ 0 ) ^ 0 , X(j)g > XfpXg 


(2.3) 



To further constrain the parameters we also require that the electroweak broken minimum 
is the deepest one at T = 0. There are several inequivalent minima whose relative ordering 
can evolve as a function of temperature. As shown in [8] this potential leads to strong hrst 
order electroweak transition over a large portion of the parameter space already at the simple 
mean field level. 

In this paper we do not consider the electroweak transition in detail, but focus on the 
properties of the DM candidate, which we take to be a singlet fermion described by the 
Lagrangian 

^DM = + sijj {gs + ifl'pTs) V’- (2-4) 

Without loss of generality, we can shift s ^ s' = s — {s) such that at the electroweak broken 
minimum {s') = 0. Since we are not assuming any discrete symmetries for the singlet scalar, 
the shift will simply reproduce all terms already present in the scalar potential with redefined 
coefficients. Hence practically we use (2.1) and (2.4), where the mass of the DM fermion at 
T < Tew is given directly by the mass parameter and the mixing between cj) and s fields 
is determined solely by In principle is a complex number, but we can always make 
it real by a chiral rotation on tp in Eq. (2.4).^ 

We use the extremization conditions to set 




= 




Diagonalizing the mass matrix 






vg.<f>s 


'^(ps 


+ Ms 


(2.5) 


( 2 . 6 ) 


we trade the parameters and with physical masses mn = 126 GeV and ms, 

^ + m%± y/ [ml - ml)2 _ 4//2^u2) , 

A = \ (jn\ + ^ y/(m^ -m|)2 -4^2^u2) . 

The scalar mass eigenstates at T = 0 are 


(2.7) 


H =/i cos/3 + s sin/3, 5 = —/i sin/3 + s cos/3, (2-8) 

where tan2/3 = —AvgL^s/{2^‘l + X^sV^ — 4A^u2). Hence, the free parameters are m^, ms, fJ-s, 
As, Xfjys, gs and gp. From equation (2.7) it follows that |/X0s| < max{mH,ms)‘^/2v. 

2.2 Constraints 

The Higgs coupling data from the LHC and Tevatron experiments constrains the mixing 
angle cos j3 and the decay width of the Higgs boson. The invisible decay width of the Higgs 
is given by 

hinv = - 2ms) + + glv^'^eimn - 2m^), (2.9) 

^Alternatively, we could rotate away either of the couplings gs,p, but then we would need to allow for 
complex Tn,jj. 
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where = 1 — 4m^/m^ and the coupling Aj/ss is given in the appendix A. We perform 
a fit to the Higgs boson signal strength data from the LHC and Tevatron experiments 
[60-66] and require that the model remains within 2a of the best fit. Moreover we check 
that the model is compatible with the electroweak precision measurements using the S and 
T parameters [67, 68]. 

We also study the perturbativity of the model as a function of the energy scale. To 
be concrete, we require that all models remain perturbative up to 10 TeV, but we find that 
many models remain perturbative up to Planck scale. The 1-loop /3 functions for the Yukawa 
couplings, gs and gp, are 

lOvrX (2.10) 

=3ffp + 5g‘^gp, (2.11) 

whereas the /3 functions of the quartic couplings read 


—24A^ -I- — 3 [Sgl + gy — 4yl) A^ 

+ g + ^gf^gy - 1 - gy) — 6yf , 

=4:X^^ ( 12 A 0 -|- 6 As) X^ps 

( 2 . 12 ) 


(2.13) 

Idn'^Pxa =13A^ -|- -|- SXsgg — 6 ( 5 ^ -|- g^) . 

(2.14) 


The (3 functions for the top-coupling yt and the weak and hypercharge gauge couplings gi 
and gy are standard, see for example ref. [69]. 

Although one can study if the model remains perturbative and maintains vacuum sta¬ 
bility all the way to the Planck scale, the question is not obviously meaningful as the approx¬ 
imation of neglecting curvature effects may not be valid [70]. Hence, we find it more natural 
to treat this dark sector as a low energy realization of some more complete theory. The 
relevant matching scale is of course completely unknown at this stage; a possible concrete 
UV completion is for example the model introduced in [71]. 

Finally, direct DM searches give an upper limit on DM scattering off nuclei. Since 
interacts with SM only through the Higgs exchange, the relevant detection channel is the 
spin-independent (SI) one. The nominal SI cross section in this model is given by: 


<^31 - 


^^NfNrn% 


TTV^ 


gj sin^/3 cos^/3 ( 


m 


H 



(2.15) 


where is the nucleon mass, /tat is the reduced mass of the WIMP-nucleon system and the 
Higgs-nucleon coupling /jv ~ 0.303. However, we do not require that V’ forms all of the ex¬ 
isting DM. We therefore define the relative DM relic density, /rei = D/i^/0.12. Consequently, 
we define an effective cross section a^ = /rei<Tgj, which is the quantity directly constrained 
by observations for a subleading DM component [13]. Currently the most restrictive upper 
limits come from the LUX and SuperCDMS experiments [72, 73]. 
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2.3 Preeze-out 

To evaluate the relative DM abundance f^ei, we concentrate to the region of parameters 
where ip is a thermal relic that decouples with the freeze-out mechanism The annihilation 
processes are ^p'lp —?• HH,SS,HS and ipip —)■ XY, where X^Y are SM fermions or gauge 
bosons. The latter are s channel processes with H or S propagators. We calculate the cross 
section in terms of the off-shell Higgs boson decay width: 

= \s^^‘^Pcos‘^f3{9svls+al) \Dh - Dsl"^ Vs Thiy/s), (2.16) 

where = 1 — s and DV = s — -|- iniaYa- We evaluate r/i(-y/s) following [14]: for 

Vs < 80 GeV and V~s > 300 GeV we use perturbative decay widths taken from [74] including 
the QGD corrections for the quark final states. In the range V^ ^ (80 GeV, 300 GeV) we use 
the tabulated results from [75], which include also four body final states. As shown in 
[76] there are large uncertainties in the range V^ G (2m7r, 2mc). In this region we use the 
perturbative spectator approach [74]. The cross section for the scalar final states is given in 
appendix A. 

Note that the annihilation cross section contains contributions both from the scalar and 
pseudoscalar couplings Qs and gp, while the cross section for DM scattering on a nucleus is 
proportional only to gs- This allows one to escape the direct detection limits by reducing gg 
while maintaining the large enough annihilation cross section by adjusting the value of gp. 

Given the cross sections, we solve the thermal relic abundance, Y = n^js (where s is 
the entropy density) from the Zeldovish-Okun-Pikelner-Lee-Weinberg equation [77, 78] 




(2.17) 

Here x = m.^/T and 




Z{x) - \j ^2 ^ ’ 

(2.18) 

where 

1/2 heS f. , T dhes\ 

^ ar ) 

(2.19) 


and ges and hes are the energy and entropy degrees of freedom, respectively. For the averaged 
cross sections, we use the integral expression [79] 


where Ki{y) are the modified Bessel functions of the second kind. Finally the relative DM 
relic density, /rei = D/i^/0.12, is given by 

Aei = 2.30 • lo9^y(xo), (2.21) 

where xq = m^/To and Tq is the photon temperature today. 

^Another possibility, requiring much smaller couplings gp^s than the ones considered here is that ip was 
ever only partially thermalized by the freeze-in mechanism. 
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3 Dark matter self-interaction 


The DM self-interaction processes are ipxj) —>■ V’V') V'V' V'V' s-iid V’V’ V’V'- The first of these 
arises from the s and the t channel, and the latter two from t and u channel diagrams with 

H and S propagators. It turns out that to obtain reasonably strong DM self-interaction S 

has to be very light, ms < 1 GeV, and the contribution of diagrams with H propagator can 
be neglected. In the interesting region the Born approximation for the self-interaction fails 
and we must calculate the cross section more accurately by solving the Schrodinger equation 
for Yukawa potential, 

v{r) = ( 31 ) 

r 

where a = cos^/I^^/dvr. In the partial wave analysis the scattering amplitude written in 
terms of phase shifts 5i is given by 

- OO 

f{0) = - j;(2/ + l)e‘'5*P,(cos0)sin,5z, (3.2) 

^ z=o 

where k = fiv = m^v/2, v is the relative velocity and Pi is the Ith. Legendre polynomial. To 
find the phase shifts we need to solve the radial Schrodinger equation, 


dr \ dr / \ r^ 


2fiV{r] 


Ri = 0, 


(3.3) 


and match the solution onto the asymptotic solution, 


lim Ri oc cos 6i ji{kr) — sind; n;(A:r), (3-4) 

r^oo 


where ji and ni are the spherical Bessel and spherical Neumann functions, respectively. In 
the numerical solution we use the method described in [29]. 

Combining all self-interaction processes the total differential cross section is finally given 
by 

^=2|/(9)-/(<r-9)l" + l/(«)iy (3.5) 

where the first term comes from the '0V'" 0V'-channels, which need to be antisymmetrized 
over the scattering angle 9, and the last term represents the V'0-channel. 

The self-interactions of DM particles are essential input for numerical simulations of DM 
halo formation. The scattering angle dependence of the differential cross section is particu¬ 
larly important when the particle trajectories are tracked through collisions. Indeed, while 
the integrated cross section is strongly enhanced in the forward region for light mediators, the 
forward collisions are unimportant as they leave the particle trajectories unaffected. Hence, 
to parametrize transport properties one uses instead the viscosity cross section 

ay = j dDsin^0^^^, (3.6) 


which in general is some function of the relative velocity of the DM particles. Moreover, we 
average the viscosity cross section over the thermal DM velocity distribution: 



d^n 

(dvrAn^) 2 




cry, 


(3.7) 
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Figure 1: Self-interaction as a function of ms for fixed a = 0.01 and Av = 10 km/sec. Solid 
and dashed lines correspond to m.^ = 400 GeV and m.,/, = 100 GeV, respectively. The shaded 
band corresponds to {ay) /m^ ~ 0.1 — 1 cm^/g. 


where Av is the radial velocity dispersion of the halo. We use Av = 10 km/sec, corresponding 
to observations of stellar motions in local dwarf galaxies [80]. In figure 1 we show the 
dependence of the self-interaction cross section on ms for fixed parameters. In order to 
be consistent with astrophysical observations, cosmological simulations require {ay) /m^ ~ 
0.1 — Icm^/g [24, 25, 27], and this region is shown as a shaded band in figure 1. 

If ms < miy/2 the LHG limits on the invisible Higgs decays constrain the parameter 
space significantly. Also the spin-independent DM-nucleon cross section, usi, which is pro¬ 
portional to l/m'g, gives a strong constraint. For a very light S the latter is by far the most 
stringent of the two. Avoiding these limits requires a very small effective coupling Qs sin /3 
when S is very light. However, Qs itself should not be small, because we want to have a large 
self-interaction rate. It follows that the mixing angle j3 has to be very small. 

In order to study the feasibility of the model, we performed Monte Carlo scan of the 
parameter space with the following priors; 

l5sl < 1, | 5 pl < 1, 0 < As < vr, -9 < logJA^^s] < 0 , 

and 

-9 < log \ii^s/ GeV] < 1.5 , -5 < log J/rs/ GeV] < 0 

1 GeV < m^ < 1000 GeV, —2 < log(m 5 / GeV) < 0. 

The resulting parameter sets, which are in agreement with the collider constraints are 
shown in figure 2 as a function of the mixing angle sin /3 and the SI detection cross section 
normalized to LUX sensitivity (left) and as a function of the DM particle mass m^p and self¬ 
interaction cross section (right). The color coding of dots specifies the mediator mass ms 
as indicated by the vertical bar to the right of the panels. To be consistent with the LUX 
limit [72] the mixing angle has to roughly satisfy a constraint 

sin/I < 10 “®(m 5 /GeV)^ , (3.10) 


(3.8) 

(3.9) 


consistent with Eq. (2.15). 
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Figure 2: Shown are the parameter sets which are in agreement with collider constraints 
and give a relic density in the range 0.8 < /rei < 1. In the right panel all points are below 
the LUX limit. Unfilled points indicate models that remain perturbative up to Planck scale. 


4 Cosmology of a light long-lived mediator 

We found that in the region interesting for DM self-interactions the mediator is light, ms < 
0.1 GeV, and its coupling to SM fermions is very small, sin/3 < 10“^. Since S decay to SM 
particles is suppressed by the small mixing angle, its lifetime tends to be very long and it 
may cause problems in the early Universe [31, 81]. In particular there will be an upper bound 
on the S'-lifetime from nucleosynthesis. 


4.1 Bound from the Big Bang nucleosynthesis 

We will first discuss the BBN bound in the case where s decays only to SM particles. Here 
S decays predominantly to the heaviest available SM fermion with a width: 


m‘i sin^/3 
^S^fJ = 87ru2 



(4.1) 


where v ~ 246 GeV and mj is the mass of the fermion. This corresponds to a lifetime 


TS 



/MeV 

V ruf 



(4.2) 


For ms ^ 2m^, S decays to muons and for 2m^ > ms > 2me the dominant channel is to 
electrons and positrons. 

When S decays to SM particles, the decay products in general thermalize quickly. An 
exception may be neutrinos, which decouple from equilibrium at T ~ 4.3 MeV (muon- 
and tau neutrinos) and T ~ 2.6 MeV (electron neutrinos) [82]. Thus, if decays take place 
around 2.3 MeV, corresponding to T5 > 0.1 sec, they may lead to an increase of the electron 
neutrino density [83] . This would increase the rate of reactions i>e + n p + e~ and u^+p 
n -|- e^, which keep protons and neutrons in thermal equilibrium and hence decrease the 
^He abundance [84]. For longer lifetimes this neutrino heating effect becomes rapidly weaker. 
Instead, for ^ 1 sec, the residual population of 5-particles starts to increase the expansion 
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Figure 3: Left: Shown are the models which are in agreement with collider constraints and 
the LUX limit and give 0.8 < f^ei < 1 as a function of the S-lifetime and the self-interaction 
cross section (uy) /m^. The vertical dotted line shows the ts = 1 sec limit and the shaded 
region corresponds to models with sufficient self-interactions. Unfilled points again indicate 
models that remain perturbative up to Planck scale. Right: The excluded regions in the 
mixing angle versus the scalar mass from BBN (blue area) and from the LUX constraints 
(yellow, orange, red) corresponding to gs = 0.4, 0.1, 0.02, respectively. The DM mass is fixed 
to 400 GeV. The dotted line indicates the expected XENONIT experiment reach [86] in the 
case gs = 0.4. 


rate of the Universe. This causes the proton-neutron equilibrium reactions to drop out of 
equilibrium earlier, which makes the ^He abundance larger. A fully quantitative analysis of 
these effects is beyond the scope of this paper. Instead we follow [85] and, to remain within 
the 2(7 limit of the observed ^He abundance, require that ts < 1 sec. 

Imposing the limit (3.10) on the mixing angle, equation (4.2) gives a rough bound on 
the 5-particle lifetime: 


TS > 0.01 


/MeVy 

\ mf J 



(4.3) 


Since for ms ^0.2 GeV the dominant decay channel is to electrons, it is clear that a light 
mediator with ms ^ 0.1 GeV, decaying to SM particles, is not consistent with the BBN 
and direct search constraints. This conclusion becomes evident also from figure 3, where 
the bounds on the mixing angle and the mediator mass following from the BBN and LUX 
constraints are shown more accurately. 


4.2 Extension with a light sterile neutrino 

An obvious solution to alleviate this problem is to introduce new states into which a light 
5 can decay to, such as a light sterile neutrino N which couples to the singlet scalar via a 
Lagrangian [59], 

(CsNN = VnSNN . (4.4) 

If rriTv < fnsl‘2^ and V mixes with the SM neutrinos with a small, but finite mixing angle sin0, 
it is possible to extend the acceptable parameter space. Let us now consider in quantitative 
detail how this mechanism works in this model. 
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First, an upper bound for the mixing sin0 is obtained from the invisible decay width of 
the Z boson, experimentally determined to be [68] 


r(Z —)• inv) 
r(Z vv) 


2.990 ± 0.007. 


(4.5) 


The invisible Z boson decay width is 


T{Z —)• inv) = r(Z —>■ Ve^e) + r(Z —)• Vt^t) + r(Z —>■ + T{Z —>■ NN). (4.6) 


If for example N mixes only with r neutrino, then 


r(z i/gt'e) = r(z = r(z w), 
r(Z —)• VrT^r) = COS^ 9 T{Z —)• Vv), 

r^\ 
ml) 


T{Z NN) = sin^ 9 T{Z ^ i^u) { 1 



(4.7) 


Now for mjv 'C mz the invisible decay with of Z is T{Z —)• inv) = 3r(Z —)> vu), as in the 
SM. This decay width is in agreement with the experimental value. The mixing modifies 
V{Z —)■ inv) only at mM ^ mz-, as shown by the red area in figure 4. 

Second, we require that the S decays to sterile neutrinos sufficiently fast and also that 
the sterile neutrinos have time to decay to SM neutrinos before they freeze-out. To be precise, 
we first require that the lifetimes of both S and N are less than 0.1 sec. The decay width for 
the process S ^ NN is the same as in the SM fermion case: 


T 


S^NN 


Vn I 1 

V = 1 - 


lOvr 


4m 


2 \ 3/2 
N 


mt 


so the lifetime of S is smaller than 0.1 sec if 


?/Ar ^ 2 X 10 


-11 


GeV\^/^ 


(—) 

\ ms J 


(4.8) 


(4.9) 


Clearly yj\f can arranged such that S decays fast enough. We have checked explicitly that 
the new annihilation channel induced by this coupling does not affect the DM abundance. 
On the other hand, the main decay channel for the sterile neutrino is to three SM neutrinos. 
The decay width for this process is 

= 292^2 "^^ ^ 

The constraint that the lifetime of N is smaller than 0.1 sec is then expressed as 

m7v(sin0)^/^ > lOMeV. (4-11) 


Hence adding a new sterile neutrino N and requiring it to decay to SM neutrinos before they 
drop out of thermal equilibrium relieves the BBN limit for ms > 20 MeV. Nominally, this 
limit excludes the region to the left from the rightmost solid line in figure 4. 

Since N's decay only to neutrinos, their late decay could lead to a significant increase in 
the neutrino energy density, potentially in conflict with the existing bounds on the amount 
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Figure 4: Shown are the excluded regions in the sterile N parameters due to the Z-decay 
width (red area) and the BBN (blue area). In the rightmost shaded region (yellow area) S 
is too heavy to provide a sufficiently strong self-interaction for 


of dark radiation [87]. We will therefore require that N energy density does not exceed that 
of one half equivalent neutrino degree of freedom at the time of their decay. Assuming that 
the freeze-out of sterile neutrinos occurred before the QCD phase transition (Tfo > Tqqd), 
we find that at T < rrie 


Pn = Tn]\}n]\} = uiN 


9es{T) 3C(3) ^3 
9es{Tio) 27r2 


O.OOSmAfT^, 


(4.12) 


which relative to the relativistic neutrinos is 


PN _ O-OOSitinT^ _ mAf 

Pu / 4 rpi T 

120 Ulf ^ 


(4.13) 


We then must require that N decays at the latest at T = mjv/10. The N decay decay width 
relative to the Hubble expansion rate is given by 




H 


10T=m^ 


0.03 


/ ruN \ 

VMeV/ 



(4.14) 


which gives 

m7v(sin0)2/^ > 3MeV. (4.15) 

Hence, the ultimate lower bound for the mass of N at given sin^ 9 comes from the requirement 
that N decays before BBN. 


5 Conclusions 

In this paper we have considered a model for self-interacting DM. The DM candidate was 
taken to be a fermion ?/) with self-interactions mediated by a scalar S. The scalar also mediates 














coupling between the hidden and visible sectors via the usual Higgs portal. We demonstrated 
that over a large region of parameter space the model can saturate the observed DM relic 
density. For the self-interactions to be sufficiently strong, the scalar mediator has to be light 
which leads to the well-known problems in the early Universe if the scalar survives until the 
era of BBN. We showed that the pseudoscalar coupling Qp makes it easier to get the correct 
DM abundance and to simultaneonsly escape the direct detection limits by reducing Qs, but 
since the DM self-interaction cross section is dominantly proportional only to gs , this does not 
alleviate the problem with DM direct detection and BBN limits in getting sufficiently strong 
DM self-interaction. We confirmed the observation [59] that this problem can be solved by 
assuming the existence of a light sterile neutrino coupling with the singlet scalar and mixing 
with the active neutrinos. The allowed region for the extended model consists of a triangle 
at 10(sin 6)-^/^ MeV < mjv < 1 GeV. 

Having established the viability of this model in the broad cosmological context, there 
are several further issues which could be addressed. A detailed investigation of various astro- 
physical probes for the DM candidate should be carried out. Also, as we have emphasized, 
the scalar sector of the model can also lead to strong first order electroweak phase transi¬ 
tion, which is a basic requirement for successful baryogenesis. We leave the more detailed 
investigation of these issues within this model for future work. 
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A Cross sections 


Here we give the formulas for the computation of the annihilation cross section for the model 
considered in Sec. 2. To make the equations more concise, it is useful to define the couplings, 
which at the electroweak broken vacuum read 


^HHH = —6nA<^ cos^ (3 — cos^ /3 sin /3 — cos (3 sin^ j3 — 2/^3 sin^ f3, 
Xhhs = QvX(j) cos^ (3 sin (3 — cos^ (3 + 2/X0s cos /3 sin^ (3 

— 2vX^s cos^ /3 sin [3 + vX^s sin^ (3 — 2/i3 sin^ j3 cos f3, 

^HSS = —QvXfp COS (3 sin^ f3 -t- 2//^^ cos^ (3 sin (3 — vX^g cos^ (3 

— ^fps sin^ (3 -|- 2uA0s sin^ /3 cos /3 — 2^3 sin (3 cos^ /3, 

^SSS = buA^ sin^ /3 — sin^ (3 cos (3 + SvX^ps sin (3 cos^ j3 — 2^3 cos^ (3, 


and 

9s,H = 9s sin/3, gp^n = 9p sin (3, gs,s = 9s cos /3, gp^s = 9p cos/3. 
Decay widths for the processes H —)• -i/iil; and S —)> V'V' 


- 2m^) 


mxR 


X 


87r 



1 - 


4m3, 




m 


+ 9v\ ^ - 


4m5, 


'ip 


‘X 


m 


X 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 
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where X = H, S and Rh = sin/3, Rs = cos/3. Cross section for the process ab —)■ cd reads 


a{ab —)• cd) = 


16TrX{s, ml, ml) 


dt \T{ab —)• cd)\^ 


where A(x, y, z) = — 2xy — 2xz — 2yz and 


t± =ml + ml± 


ml + 


' X{s,ml,ml)X{s,ml,ml) 


X{s,ml,ml 


ml + 


X(s,ml,ml 


The squared amplitude, averaged over the initial states and summed over the final states, 
for the DM, annihilating to two scalars, hi,hj = H,S, or equivalently for hi and hj 
annihilating to -pair is given by 


= (l-f 


fji{ss) ^ ^ ^ t2jn{su) ^ 2^fjn{tu) 

IJ IJ IJ IJ IJ IJ 


where 


dss) _ , 


= {2-dij) {s-4ml) 


\k=H,S 


Xijk9s,k 

s-ml + imfcTfc 


(A.IO) 


, mjkdp.k 1 

+■5 / -> 

= (2 (m^ -ml{t + u)- mfm] + tu) 

+2 (m^ + mliu- 3t) - mlfm] + tu) (5's,i5'pj + 5'p,i5s,i) 

-2 (7m^ + - 4mj +9t + u) + m^mj - tu) gp^ig^j 

T8?tt. sgs^igs,jgp,igp,j) {t ^ip) i 

T(uu) ^ ^2 (^4 _ ^2 _ ^2^2 ^2 ,^2^ 

+2 (m^ + ml{t- 3u) - mlfm] + tu) {g^igl^j + glpdl^j) 

-2 (7m^ + m^(-4mi - 4mj + 9ii -|- t) + - tu) g^pg^j 

T8?tt. sgspgs,j9p,igp,j) '^-ip) > 

T(^t) ^ _ Yh _ 

X [[2mp,{t - u)gspgs,j + 2m^(4m^ - 2ml - 2m| + 3t + u)gppgpj) gs,k 
+ {2mpj{m‘f -ruj + s)gspgp,j + 2mp,{-ml + m] + s)gppgs,j) gp,k) , 
xi^t) ^ _ 

X {{2mp,{u - t)gspgs,j + 2m^(4m^ - 2ml - 2m| + 3u + t)gppgpj) gs,k 
+ {2mpj{m‘f -ruj + s)gs,igp,j + 2mp,{-ml + m] + s)gppgs,j) gp,k) , 


(A.ll) 


(A.12) 


(A.13) 


Tif= E 


(A.14) 


13 - 



(A.15) 


= (2 (m^ -w?^{t + u)- mfm] + tu) gl^gl^ 

-2 (m^ + m\{t- 3m) - + tu) glig^j 

-2 (m^ + m^(M - 3t) - + tu) gl^ig^j 

-2 (9m^ + mli-Anii - Am] + 3{t + m)) - mjm] + tu) gl^tg^j 
+3m^sgs,igs,jgp,igp,j) {t - m^)“^(M - m^)"^ 
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